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1.  Introduction 

Multivariate  extreme  value  distributions  have  been 
studied  by  many  authors,  and  their  contributions  are 
summarized  by  Galambos  [1]  and  Resnick  [2].  The  pur- 
pose of  this  paper  is  to  obtain  some  necessary  and  suffi- 
cient conditions  for  domains  of  attraction  of  the  multi- 
variate extreme  value  distributions.  The  joint  asymptotic 
distribution  of  multivariate  extreme  statistics  is  also  ob- 
tained. To  study  multivariate  extreme  value  distributions 
and  their  domains  of  attraction,  Sibuya  [3]  introduces 
the  notion  of  a dependence  function  which  is  also  used 
by  Galambos  [1].  A dependence  function  or  copula  is  a 
useful  notion  to  construct  a family  of  joint  distributions. 

In  this  paper,  basic  arithmetical  operations  are  always 
meant  componentwise  (see  Galambos  [1],  Chapt.  5). 

Let  (Xi;,X2;,...,Xij),y=l,2,...,rt,  be  a sample  of  size  n, 
of  a /: -dimensional  random  vector  with  a distribution 
function  F{x).  The  /-dimensional  distribution  function 
of  the  components  Xj^,Xj^,...,Xj.  will  be  denoted 

We  shall  also  use  the 

notation  Fj,..j.(T;,,...,T;,)=Fy(o(X2(,))=P(X,,>T2,,...,X2>T;,). 
For  /:=1  and  pE(0,l),  let  F~'{p)='\ni{x:F{x)^p}. 

Let  Z„=iZu,...,Zk„),  where  Z,„=max{X,i,...,X,„ }, 
i=\,2,...,k,  and  let  us  call  Z„  a multivariate  extreme 
statistic. 


If  there  exist  a„>0,b„E.R\n=\,2,...{a„>0  means 
a,„>0,/=l,—,^)  such  that  {Z„—b„)la„  converges  in  distri- 
bution to  a random  vector  U with  a nondegenerate  distri- 
bution H (i.e.,  all  univariate  marginals  of  H are  nonde- 
generate), then  F is  said  to  be  in  the  domain  of  attraction 
of  H , FGD{H)  by  symbol,  and  H is  said  to  be  a multi- 
variate extreme  value  distribution.  The  convergence  in 
distribution  is  equivalent  to  the  condition 

\\m  F”ia„x+b„}=H (x)  (1) 

/?— »» 

for  all  X,  because  multivariate  extreme  value  distribu- 
tions are  continuous. 

We  shall  need  the  following  lemma  to  prove  a propo- 
sition in  Sec.  2. 

Lemma  1.1  Equation  (I),  is  equivalent  to 

lim  n [1  — F(a,;t:-i-fe„)]=-log  H{x) 

for  all  X such  that  0<H  (jc)<L  (See  Marshall  and  Olkin 
[4].) 
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2.  Domains  of  Attraction 

For  any  ^-dimensional  distribution  F , 
D,(y)=F{FV(y^) Fl'iyk)),  y=(y„...,y*)e(0,l)* 


^HCv;/')=yi-yr  exp{X(-iy  S 

i-2  \sj  ^<-<jiSk 

and  DH{y\r)  is  a dependence  function  of  a multivariate 
extreme  value  distribution. 


is  called  the  dependence  function  of  F.  In  this  section.  Proof.  It  is  easily  seen  that  for  all  5>0, 
we  derive  necessary  and  sufficient  conditions  for  do- 
mains of  attraction  in  terms  of  the  dependence  function.  5^^y(/)(Cy7(/))'^0=^A/)(yA/))- 


Proposition  2.1  Let  F be  a k-dimensional  distribution 
and  let  H be  a multivariate  extreme  value  distribution 
with  univariate  marginals  //,,  i=  1 ,...  ,k.  Then  the  follow- 
ing statements  are  equivalent: 

1}  F^DiH). 

2)  F.GDiHi),  i=l,...,k,  and 


From  Theorems  5.3.1  and  5.2.4  of  Galambos  [1],  we 
have  the  result.  □ 

Example  2.1  (See  Examples  5.2.2  and  5.2.3  of  Galam- 
bos [1].)  For  a Mardia’s  distribution 

F(jr,,jC2)=l  -e“^'  -e-Me"'+e"2- 1)'', 


V\m  s[l-DF(y''')]=-\og  Dniy)  for  all  yG(0,l)*. 
3)  F,E^D{Hi),  i=\,...,k,  and 


lim 

xTi 


l-D^Cv'-Q 

. 1-^ 


= -\og  Dniy)  for  all  yG(0,l)\ 


4)  F,^D{Hi),  i=\,...,k,  and 

Proof.  The  proof  is  straightforward  from  Lemma  1.1, 
Theorem  5.2.3  and  Lemma  5.4.1  of  Galambos  [1].  □ 

Proposition  2.2  Let  F be  a k-dimensional  distribution 
and  let  H be  a multivariate  extreme  value  distribution 
with  univariate  marginals  Hi,  i=\,...,k. 

(A)  FGD{H)  if  and  only  if  FiE.D{Hi),  i=  \,...,k,  and 
the  functions 


for  each  fixed  vector  J{i){i>\)  and  for  all  are 

finite,  and  the  function 


k 

£>wC)’T)=yi-y7  exp{2(-l)'  2 dmiym)} 

i-2  is;, 


is  a dependence  function  of  H. 

(B)  The  following  inequalities  hold. 

D/,  (y  ;2r+ 1 )<D„  (y  )<Dh  (j  ;2r ), 

for  a nonnegative  integer  r,  where 


T>fCVl,>’2)=>'l+>’2-l  + 


and 


nDF{yl'\yf'')=n 


1 -y,  1 -y2 


1 1 
-I- 


Ll-y,'"'"  l-y2‘'” 

^_(logy,)(logy,)^ 
log  yi-i-log  y2 

Thus,  by  Proposition  2.2  we  have  F G D{H),  where 
T>H(yi,y2)=yry2  exp 


(logyi)(logy2) 
logyi-t-logy2  J’ 


H{x\,  X2)=A  (X|  )A  (jC2)exp{  1 /(e''-t-e^^) }, 
and  yl(jr)=exp(— e”0- 

Proposition  2,3  Let  F and  G be  k-dimensional  distri- 
butions and  let  H be  a multivariate  extreme  value  distri- 
bution. 

1)  IfFGEDiH),  then 

2)  If  F^D{H),Gi^D{Hi),i=\,...,k,  and 


lim  I where  l={\,...,\), 

yTi  l-Doiy) 


then  GGD{H). 
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3.  Marginally  Independent  or  Perfect 
Dependent  Multivariate  Extreme  Value 
Distributions 

Let  //  be  a multivariate  extreme  value  distribution 
with  univariate  marginals  //,,  Let 

and  //*(x)=min{//,(jr,),  }, 

then  it  holds 

H,{x)<H{x)<H\x) 

for  all  xG/?*.  Both  bounds,  //.  and  H*,  are  multivariate 
extreme  value  distributions.  Characterizations  of  these 
distributions  are  obtained  by  Takahashi  [5], 

In  the  bivariate  case  Sibuya  [3]  obtains  necessary  and 
sufficient  conditions  for  FGDiH*)  and  FE:D{H*).  In 
this  section  we  generalize  his  results. 

Proposition  3.1  Let  F be  a k-dimensional  distribution 
and  let  Hi  be  a univariate  extreme  value  distribution, 
/=!,... Then  the  following  statements  are  equivalent: 

1)  F&D{H.). 

2)  FiE.D{Hi),  i=\,...,k,  and  there  exists  yG{0,\Y  such 
that 

\\m{Dr{y'"')T=yiyr-y,. 


3)  FiGDiHi),  i=\,...,k,  and 


lim- 

>Tl 


l-£>F(yl) 

l-y 


=k. 


By  Proposition  3.1  4)  we  have  FElD{H»),  where 


Proposition  3.2  Let  F be  a k-dimensional  distribution 
and  let  Hi  be  a univariate  extreme  value  distribution, 
i=\,...,k.  Then  the  following  statements  are  equivalent: 

1)  FGDiH*). 

2)  FiGDiHi),  i=\,...,k,  and  there  exists  yE(0,l ) such 
that 

lim(Df(y'^T)r=y. 


3}  FiEDiHi),  i=\,...,k,  and 


y-t\  l-y 


Proof.  The  proof  is  straightforward  from  Theorem  3.1 
and  Corollary  3.1  of  Takahashi  [6].  □ 


4.  Joint  Asymptotic  Distribution  of  tbe 
Multivariate  Extreme  Statistics 

In  this  section,  we  show  the  joint  asymptotic  distribu- 
tion of  several  multivariate  extreme  statistics  along  the 
arguments  in  Sec.  2.3  of  Leadbetter  et  al.  [7].  For  sim- 
plicity we  shall  consider  the  bivariate  case. 

Let  (Xi,T|),...,(X„,y„)  be  a sequence  of  independent 
random  vectors  with  common  distribution  F.  The  order 
statistics  of  the  components  will  be  denoted  by 


4)  FiE.D{Hi),  i=l,...,k,  and 


; and 


lim 

yTi 


1-T>r(yl) 

l-y* 


= 1. 


For  /=0,l,...,r— 1,  define 


Zn-i—{X„-ij, , Yn-i-.n  ) 


Proof.  The  proof  is  straightforward  from  Theorems  2.2 
and  4.1  and  Corollary  2.4  of  Takahashi  [6].  □ 


and  let  us  call  Z„  , an  (/-i-l)-th  multivariate  extreme 
statistic. 


Remark.  If^=2,  we  have  the  same  result  as  Proposition 
3.1  by  Corollary  2.2  of  Takahashi  [6]. 

Example  3.1  (See  Example  5.2.3  of  Galambos  [ 1 ].)  For 
the  Morgenstern  distribution 

F(xi,X2)=1  — e — e“^^-i-e“"''“'"41-i-a(l  — e"^‘)(l 

(yi,>’2)=yry2[l+a(l  -yi)(l  -y2)] 

where  — l<a<l,  and 


Proposition  4.1  Suppose  that 

P{(Z-b„)la„^x}^Hix) 

for  some  nondegenerate  distribution  H.  Then,  for 

Xi={Xi,yi)>X2={x2,y2), 

P{{Z„- b„ )la„<Xu  {Zn-\- bn fUn^Xj }^Hj (X| ,X2 ) 

where 


lim- 

>-Ti 


l-D/r(y,y) 
1 _„2 


Hj{xuX2)=H{x2){\-v\og  H (xi)-\og  H {X2) 

+[\ogHdxi)-logHi{x2)+ih{Xi,y2)-h(x2))] 

X [log  //2(y  I ) - log  H2(y2)+{h  (x2,y  1 ) - /?  (X2))] } 
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and  hix)=\im„^a=  F{a„x+b„). 

Proof.  Define 

5o=#{/  \X->ai„Xt+bi„  or  F^a2n>’i+^2n  J=h2,...,n.}, 
S['=#{/  \ ainX2+b\„<Xj<ai„Xi+bx„  and 
Yj^a2„y2+b2n  J=  1 ,2,...,n . }, 

52=#{/  \Xj<a\nX2+b]„  and 

a2ny2+b2r,<Yj<a2„yx+b2„  J=  1 ,2,...,n . }, 

Sx2=#{j  I a„X2+b„<{Xj,Yj)<a„Xx+b„,j=l,2,-,n.} 

then,  we  have 

P bfi')^afi—X\j  {Zn~\  bti')lan—X2y 
=P  {5o"=0,  5, "2=0,  1 , 1 } 

+F{5o"=0,  SiVl,  5,'’=0.  52"=0}. 

On  the  other  hand,  by  using  Theorem  5.3.1  of  Galambos 
[1],  we  can  evaluate  the  asymptotic  probabilities  of  the 
evens 

{5u"=/,  S^2=j,  S^=k,  S^=m  } 
for  i,j,k,m=0,\.  Thus  we  have  the  result.  □ 

Corollary  4.1  Suppose  that 

P {{Z„—b„  )la„<x  }A//(x) 

for  some  nondegenerate  distribution  H.  Then,  for  fixed 
rS:l  and  Xi>-">Xr 
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\P{{Z„-b„)la„<Xx,...,{Z„-r^x-bn)la„<Xr) 

-P  {(Z*-/3„)/a„<x, (Z*_,+|  )la„<Xr } I ^0, 

as  n-^^, 

where  Z*_,  is  the  {i-\-\)-th  multivariate  extreme  statistic 
from  the  distribution  H,  i=0,...,  r— 1,  and 
H''{a„x+fin)=H {x),  n=l,2,... 


Example  4.1  Let  F be  the  bivariate  normal  distribution 
with  the  correlation  coefficient  less  than  one.  Then 

\ P {{Z„-b„)!a„<Xu...,{Z„-r+x-bn)la„<Xr} 

-P{(Z;-(log  «)l<Xi,...,(Z„‘_,.„-(log  r)1<x,  } I ->0, 
as  «— >00, 

where  Z*_,  is  the  (t+l)-th  multivariate  extreme  statistic 
from  the  bivariate  exponential  distribution  whose  mar- 
ginals are  equal  to  the  standard  exponential  distribution 
and  they  are  independent.  For  the  univariate  case,  it  is  a 
well  known  result  (see  Weissman  [8],  Theorem  3). 
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